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Abstract 

The finite temperature Casimir effect for a scalar field in the bulk 
region of the two Randall-Sundrum models, RSI and RSII, is studied. 
We calculate the Casimir energy and the Casimir force for two par- 
allel plates with separation a on the visible brane in the RSI model. 
High-temperature and low-temperature cases are covered. Attractive- 
ness versus repulsiveness of the temperature correction to the force is 
discussed in the typical special cases of Dirichlet-Dirichlet, Neumann- 
Neumann, and Dirichlet-Neumann boundary conditions at low tem- 
perature. The Abel-Plana summation formula is made use of, as this 
turns out to be most convenient. Some comments are made on the 
related contemporary literature. 

PACS numbers: 03.70. +k, ll.lO.Kk, ll.lO.Wx, 04.62.+V 



1 Introduction 

Inspired by the Randall-Sundrum models [Ij, there has recently been con- 
siderable interest in the Casimir effect in higher-dimensional space. We may 
recall the characteristic features of this model: in the first variant, called 
RSI, one assumes that we are living on a (3+l)-dimensional subspace called 
a 3-brane, separated from an additional hidden brane by a bulk region. Only 
gravity is assumed to propagate in the bulk. The extra dimension is a circle 
with radius r^, represented by a coordinate in the range — vr < < vr. 
The hidden and the visible branes are located at = and = vr respec- 
tively. Imposition of Z2 symmetry means that the points {x^, 0) and {x^, —0) 
are identified. In the second variant of the model, RSII, the hidden brane 
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is sent to infinity. The major difference between the RS model and other 
higher- dimensional models lies in the warp factor e~^'^'''=l''^l in the metric 

ds^ = e-^'^'^l'^lr^^^rfx^rfx'^ - rld(j)^, (1) 

where 77^1^ =diag(l,-l,-l,-l) is the Minkowski metric of flat spacetime and 
is a constant of order Mpi, the Planck mass. The warp factor plays an 
important role helping to solve the hierarchy problem without introducing 
additional hierarchies. 

Whereas in the original RS model the fields of the standard model (SM) as 
mentioned were confined to the visible brane only, the possibility of having 
additional fields in the bulk was soon investigated, beginning with scalar 
fields [21 [3] introduced to stabilize the inter brane distance. Subsequently, 
the possibility of having other fields, such as fermion field [H [5], gauge fields 
[ni [71 E], or even the full assembly of SM fields [S], was investigated. Tests 
of Newton's law at short distances may show aberrations at short distances 
(cf., for instance, Refs. [9]). 

Only recently have there appeared papers on the Casimir effect in the 
Randall-Sundrum models. To our knowledge the first group working on this 
was Frank et ai, publishing two papers fiUi [TT] on the Casimir force in 
both the RSI/RSII and the RSI-q/RSII-q models. Here RSI-q and RSII-q 
refer to generalizing the 3-branes on RSI/RSII to (3 + g)-branes (this kind of 
generalization will however not be dealt with in the present paper). There are 
also two papers of the group of Morales- Tecotl et al. [12], focusing on RSI- 
q/RSII-q. While Frank et al. used zeta function regularization to calculate 
the Casimir force. Morales- Tecotl et al. used a Green's function formalism. 
A delicate point is that their results are seemingly in conflict except in the 
Minkowskian case. We shall comment on this point later. 

A third class of papers are those of Cheng [131 [IH US]- These papers, 
except from Ref. [15], assume zero temperature. Moreover, Elizalde et al. 
have recently studied repulsive Casimir effects from extra dimensions for a 
massive scalar field with general curvature parameter pi3j. 

Among previous papers in this research field, the ones most similar in 
nature to that of the present work are those of Teo [171 [HI [TH [20l [21] . She 
calculated the temperature Casimir force; not, however, the Casimir free 
energy. 

As already mentioned we will assume a scalar field $ in the bulk. Both 
the RSI and the RSII models will be considered. Formally, the expressions 
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pertaining to the RSII case may be derived by letting ^ oo in the RSI 
expressions. The use of a scalar field makes of course the situation more 
unphysical that what would be the case by assuming an electromagnetic 
field in the bulk. But we avoid the complications arising from photon spin in 
higher- dimensional spacetimes (for some recent papers in that direction see, 
for instance, Refs. [21 [231 El E51 EH EZ] ) ■ The confiict is not even resolved 
when we have only one extra spatial dimension and spacetime is flat. On one 
hand we have e.g. Poppenhaeger et al. [27] and Pascoal et al. [26]. They find 
the electromagnetic Casimir force by multiplying the scalar field expressions 
by a factor p (to account for the possible polarisations of the photon) and 
subtract the mode polarised in the direction of the brane. On the other hand 
we have Edery and Marachevsky [25] who start out with decomposition of 
the five-dimensional Maxwell action. This confiict in not the central issue of 
this paper, and we avoid it by consider scalar fields only. 

Our main purpose will be to calculate the Casimir free energy and the 
Casimir force for the RSI model, when there are two parallel plates with 
separation a on the visible brane. This is the piston model. Our main focus 
will be on the following points: 

1) The calculation is given for arbitrary temperature T, and the low- 
temperature and high-temperature limits are thereafter considered. The 
attractiveness versus the repulsiveness of the temperature corrections for 
different boundary corrections at low temperature are of definite physical 
interest and are therefore pointed out. We regularize infinite expressions by 
using zeta functions and the Abel-Plana summation formula, as this formula 
turns out to be better suited to the problem than the more commonly used 
Euler-Maclaurin formula. 

2) We assume Robin boundary conditions on the physical plates at a; = 
and X = a. Usually, one has been considering the more simple Dirichlet 
conditions when working on this kind of problems, although very recently 
the Robin considerations have begun to attract attention [281 121] • 

As an introductory step, we consider in the next section the partition 
function and the free energy of a bulk scalar field. We discuss the distinc- 
tion between even and odd fields, and consider also the mode localization 
problem. After a brief survey of the Abel-Plana formula in Section 3 we con- 
sider in Section 4 the Dirichlet-Dirichlet (DD), Neumann- Neumann (NN), 
and Dirichlet-Neumann (DN) boundary conditions in fiat space, at finite 
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temperature without extra dimensions. Our main topic, the temperature 
RSI case, is covered in Section 5, where the Casimir free energy and force are 
calculated for the different boundary conditions. A brief treatment of the 
RSII case is given in Section 6. 

It should be recognized that the warp factor is an important element 
in the present problem. One might analyze instead the analogous higher- 
dimensional cases taking spacetime to be flat. Considerable interest has 
been devoted to this simpler variant of higher-dimensional Casimir theories 
in recent years. See, for instance, Refs. [HI [291 [25l [301 [311 [321 [26j , cind further 
references therein. 



2 Free energy of a bulk scalar field 

To find the partition function for a non-minimally coupled scalar field $ with 
mass m in the RSI model, we follow a Kaluza-Klein reduction approach [33] . 
starting from the Lagrangian density 

C = Q^M^a^'^ - + Ci? + Chid5(^) + c^-Jiz - j . (2) 

Here G'=det Gmn (with M,N = 0,1,2,3,5) is the determinant of the 5D 
metric, R is the 5D Ricci scalar, ( is the conformal coupling, and Chid/vis 
are the boundary mass terms on the branes. Throughout the article we use 
h = c = ks = 1- We have introduced above a new position coordinate z such 
that \z\ = (e'^''"'^'^' — l)/k, implying that Zr = (e'^^'^" — l)/k. It is convenient 
to introduce also the quantity A{z) = 1/(1 + k\z\). The partition function 



Z = j exp j d^xdzC^ 



(3) 



can now be calculated, making use of the Euclideanization = x*, x° = r 
ix^ . A partial integration yields 



j exp - j (fMz]^^A^{z) (f + $ 



(4) 



Here 

p2 = f,^^-^^^,, (5) 
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where t]'^'^ = —S'^'^ is the metric in the coordinates x^, and 

= A-''{z) [-d,A^{z)d, + A'{z) (m^ + CR + ChuS{z) + c,,,5(^ - z^))] . 

(6) 

We now expand $ in the eigenf unctions Xp{^^) ipNiz) of and M| 
respectively, 

'^'l^. ^) = cjv(p)Xp(5)^Jv(2;), (7) 

N,p 

where (r = it) 

p\p{x) = -(5^ + 9^ + dl + dDxpix) = p\,{x), (8) 

M^ijNiz) = MlijM{z). (9) 
The eigenfunctions are normahzed as 

d^xXp{.S:)xp'{,x) = 5pp>, (10) 

dz'4)N{z)A^{z)i)N,{z) = 5nN'- (11) 
The partition function now takes the form (an unimportant factor omitted) 

Z=X{ {MI+pY""^ (12) 

A/jv,p 

where the sum goes over all eigenvalues of and p. Our next step is now 
to identify M^r and p. 

2.1 Eigenfunctions and eigenvalues for fp' and 

We start from Eq. ([H]), assuming Robin boundary conditions on the physical 
walls, 

(l + /3o9,)Xp(£)|,=o = 0, (13) 

(l-/?,9.)xp(x)U, = 0, (14) 

with constants (3^ and (3a referring to a; = and x = a. The forms above are 
as in Ref. ^28j. Dirichlet and Neumann boundary conditions correspond to 
/5 = and (3 = oo, respectively. We assume eigenfunctions of the form 

Xp{x) = Ne'^''^+''yy+''^'^ cos{Kx + a), (15) 
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with eigenvalues 



p^ = e^ + kl + kl + kl (16) 



For a bosonic field at temperature T the Matsubara frequencies are 

Q = 27rT/, / G Z (17) 
Equation f|T3|) leads to the following constraint on a, 

cosq; = — :, (18) 

whereas Eq. f|T4l) yields F^ikx) = 0, where 

i^x(fcx) = (1 - klPoPa) sm{k^a) - k^{j3o + (3a) cos{k^a). (19) 

Consider next the eigenf unctions ilj]sr{z). We insert the expression into 
Eq. iQ, take into account that the Ricci scalar for the RS metric is 

R=-20k'^ + 16k{S{z)-5{z-Zr)), (20) 

and change the position coordinate in the bulk back to y using d/dz = 
A{y)d/dy. Then, 

^J^(y) - Ak^p'j^iy) + {Mf^e^'y - (m' - 20Ck')) Mv) = 0- (21) 

The solution is (we consider the region < y < nrc only) 

My) = ^ U (^e^A + K{M,)K (^e^')) , (22) 



where ly = ^^4 + (m/k)'^ — 20( and Cat a normalization constant. This is 
the same result as in Ref. except that we include curvature {( ^ 0) in 
our model. 

One should now distinguish between even fields satisfying ip]y{—y) = 
tpN{y), and odd fields satisfying ipiy{—y) = —ipiy{y). The behavior may be 
summarized as follows: 

• Even scalar fields obey the Robin BC on the branes. If the field is mini- 
mally coupled (C = 0) and there is no mass boundary term {cbrane = 0), 
the boundary condition reduces to the Neumann BC, ip'i^{y)\brane = 0. 
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• Odd scalar fields obey the Diriclilet BC on the branes. 
These two cases may be combined: we introduce the two functions 

J^-(^) = (2 - (fc/?b.anc)-^) J.(^) + ZJI{Z) 
^brane(^) ^ (2 - (A;/5b,anc)-')nW + ^n'W, 

and let now z mean z = e^'^^'^M^ / k (not to be mixed up with the coordinate z 
in Sect. 2), and d = e~^'^'^''. Then we can write the general BC as F]s[{z) = 0, 
where 

F^iz) = 3f\zd)yf{z) - 3T{^)y':\^d). (24) 

This is in accordance with Ref. [3S] in the case of minimal coupling, if we 
choose Chid = -Cvis = 2a /k. 

Special attention ought to be given to the massless case, = 0. For 
fields with — 20(k'^ ^ there is no solution of Eq. fl2Tl) with Mn = 
satisfying the Robin BC on both branes. For an even field with — 20(^k'^ ^ 
with no boundary mass term the situation is different, as t/'o =const is a 
solution of Eq. fl2T|) and also satisfies the boundary condition which in that 
case is the Neumann BC. The Mat = case has important consequences for 
the Casimir force from a bulk scalar field. This is related to the localization 
problem for the Kaluza-Klein modes in general. In RSI, the massless mode 
is localized near the hidden brane at ?/ = 0. In RSII the situation is reversed, 
as the massless mode is localized near the visible brane at |/ = and the 
massive modes are delocalized. The reader may consult Refs. |36] for a 
discussion as to what weight to be given to the massless modes in RSI due 
to the fact that it is localized near the hidden brane only. 

2.2 Approximate expressions for the masses 

We assume d = e~'''^^'' -C 1 but keep z = e^'^'^'^M^ /k arbitrary, to find 
convenient approximative expressions for the Kaluza-Klein masses. As in 
this case i^"''^^{z) <^ yl^"'^'^ (z) , the equation Fjq{z) = reduces to 

jriz) = 0. (25) 

The situation can be divided into two classes: 

(i) For Dirichlet BC {Pbrane = 0) it follows that we need the zeros of J2{z). 
Making use of the large-z approximation Ju{z) ~ (2/7r)^/^ cos{z — ^utv — W) 
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we find that the expression 



M;v = kTie-"-'^ (^N+^u-^j, N=l,2,... (26) 

is useful for practical purposes. 

(ii) For non-Dirichlet BC {(3hrane 7^ 0) we obtain from Eq. fl23|) 

j,{z) = (2 + z/ - {kl3Y^)J,{z) - zJ,+,{z), (27) 

leading approximately to 

Mm = kne-'-"^^ (^N + ^u-^y N=l,2,... (28) 

The formula is good for kjSbrane > 1 and becomes better for higher A^. As 
an example, choosing kjSbrane = 10^, d = 10~^^, z/ = 2, we find the numerical 
error of the zeros to be about 4% when = 3 and around 1% when N = 5. 
As we will see later the first (i.e. smallest) values of Mn are the most 
significant for the Casimir force in RSI. 



2.3 Two expressions for the free energy 

From Eq. (fT2|) we obtain for the free energy 



F = -T\nZ = -TV^ ^ 



r-ljiv^T. E HMl + eUkl + kll (29) 

where k']_ = k'^ + k"^, V± is the transverse volume, ei = 2ttTI, and the summa- 
tions over kx and Mjy go over all real zeros of the functions Fx{kx) (Eq. (|T9l) ) 
and Fn{z) (Eq. ([24])). 

By making use of the zeta function 

C(-)= E ^ E {Ml + eUkl + kl)-\ (30) 

l=-oo ^ ' Mti,k^ 

following [3lj, we can re express F as 

^ \^l:^''as) „, (31) 



2 ds 



s=0 



where n is an arbitrary parameter with dimension mass. 

We now derive the classical expression for F using that the Mellin trans- 
form of b~^r{z) is e~^^, i.e., 

1 

f-^e-'^'dt. (32) 







Applying the Poisson summation formula (details omitted) we can then de- 
rive 

40F (27r)^£t .t-^Jo ^ ' 



Further manipulations lead us to the desired expression 

^ = ^^Vl^ 5: In hsinh (^^Mj, + kl + kl 
One may note here that a boson with energy Ep contributes with (/? = 1/T) 



(34) 



oo 

"'h -2sinh(/3E,/2) ^^'^ 

to the total partition function [35]. Summing over all energies we obtain 
the classical expression corresponding to Eq. flMl) . The expression of the 
free energy of a scalar bulk field is equal to that of bosons with energy 
E^ = Ml + kl + k\, where Mn is the masses found in Section 12.11 For 
Mn = this is free energy of a scalar field in Minkowski (i.e. flat) spacetime 
without extra spatial dimensions. By letting T — we find the zero-point 
energy 

^ = ^^1^ E y^l + kl + kl. (36) 

Again, we observe that Mn = 0-term in the sum correspond to the familiar 
expression for 3-1-1- dimensional Minkowski spacetime. 

Another expression for F can be derived which in our context is more 
useful, in view of our application of the Abel-Plana summation formula later 
on. We start from the expression fl3Tl) . introduce a generalized polar coordi- 
nate transformation along the same lines as in Ref. and integrate over 
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all angles. We then obtain (the limit s ^ is understood) 

TVx_ d 



Arc ds 



l=-oo MmMx 



(37) 



where C is defined as 

C = M'^ + e^ + kl. (38) 

The integral is solved using the variable change x = r'^/C and leads essen- 
tially to the Beta function B{q,v) = r(g)r(t>)/r(t> + q). We obtain 



TV± 



r(-l) E Y.^^N + e^ + kl 

l=—oo Misr,kx 



(39) 



This is the finite-temperature form that we will use below. The corresponding 
zero-temperature form is found by a limiting procedure to be 



E 



(40) 



From now on we will set = 1. Hence F and P (force) refer to respec- 
tively energy, free energy and force per unit area of the physical plates. 

2.4 The piston model 

Before finding explicit expressions and specifying BCs we introduce the piston 
model. The model has attained a great deal of attention [3HI EHl [29l HOj HH 
|2T] . We introduce the piston (Figure [I]) with the same notation as in Chapter 
4.3 of [l2]. Instead of only using the free energy Fi of cavity / as the Casimir 
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Figure 1: Illustration of the four cavities in the piston model. 

free energy, we use 

Fpiston = F\d) + F'\X -a)- F''\X/r,) - F^^(X(1 - l/r,)). (41) 
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Initially the system is in an unstressed situation where the cavities have size 
X/t] and X{1 — 1/ rf). Then we shift the middle plate so that the lengths of the 
two cavities are a and X — a; the system is now in a stressed situation. The 
Casimir free energy is the sum the free energies of two cavities in the stressed 
case (I and II) minus the free energies of the cavities in the unstressed case 
(III and IV). The constant is ~ 2, characterizing the unstressed situation. 
In the end we let X ^ oo and effectively remove the rightmost plate from 
the setup. In the piston model all terms independent or linear in a vanish, 
hence from now on we will discard all such terms. 



3 Casimir free energy and force: Initial re- 
marks 

3.1 The Abel-Plana formula 

We want to find a more explicit expression for the Casimir free energy, one 
we can evaluate numerically. Thus all the summations over and Mn from 
Eq. ( l39l) need to be taken care of. Instead of using Eq. fl39l) we look at the 
complex function 

m = --r{s) J2 J2(^N + el + kl)-% (42) 

/=— oo MM,kx 

which reduces to the free energy in Eq. (1391) when s = — 1. The function 
F{s) is well-defined for large, positive Re(s) and we analytically continue it 
to the whole complex plane. Together with F{s) we will use a variant of the 
Abel-Plana formula \M EH] 



y vr/(^„) /(O) 

4-f 1 + (l/^n) sinz„cos(z„ + 2a) 2 1 - /3o/a - /3a/a 

n-i s / 

(43) 



oo roo 



" . >^U (A,/a+l)(/3a/a+l)'' 



especially suited for plates with Robin BC. Here, Zn denotes the n'th zero in 
the right half of the complex plane of the complex function Fx{z = akx) in 
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Eq. fll9p . From Eqs. fllSp and 0191) we can find the relation 



1 + (l/z„) sinz„cos(^„ + 2«) = 1 - y ^ - , (44) 
By choosing 

f(') = Z « + + ^Va^) f 1 - E I ' (45) 

V j=0,a 



the left hand side of Eq. fj43l) matches the sum over in -F(s). The notation 
Y2j=o a i^sans there are contributions from both the left (j = 0) and the right 
(j = C') plate. 



3.2 Application of the Abel-Plana formula 

We can divide the free energy F at arbitrary temperature T into two separate 
parts, F = F{Mn = 0) + F{Mn > 0). For a massive scalar field there is 
no massless mode (Mjv = 0) at all. For a massless field, even and minimally 
coupled, there is an Mn = mode. Recall that F{Mn = 0) yields the same 
expression as the free energy of the massless scalar in Minkowski spacetime. 
To find the Casimir energy and force for such a field one can simply add the 
massless mode term. It is natural therefore to analyze the Mn = mode 
separately. The formal expressions are divergent, and will be regularized by 
the use of zeta functions. 

Let now be arbitrary. Insert the expression (H5|) into P3|l . and divide 
the sum into three separate parts as indicated by the underlines 1,2, and 3. 
We do not give the details here, as the formalism is analogous to that of 
Ref. [2B], pertaining to the zero temperature case. The free energy can be 
written as the sum of three parts: one part F^p as the contribution when no 
plates are present, one part Fj as the vacuum free energy along the transverse 
directions induced by the plates at = and Xa = a respectively, and a 
remaining part AF. Thus 

F = Fnp 

j=0,a 

The two first terms do not refer to the gap width a, or are linearly dependent 
on a, and do not contribute to the free energy in the piston model. The last 
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term AF, henceforth called simply F, is the term of physical importance. It 
is precisely the term corresponding to underline 3 in Eq. We give this 

expression explicitly: 

T 



[zVa^-(M^ + 6n]- 

J a fWT? (/3o/a-l)(/3a/a-l) 2. _ ^ 

-oo Mjv "''^'V^^iv+^i (/3o/a+l){/3a/a+l)'^ 

1 - y . (47) 

We can now use this expression as basis for discussing special cases: DD, NN 
and DN boundary conditions. We first consider flat space with no additional 
spatial dimensions. 



4 DD, NN and DN boundary conditions in 
flat space with no extra dimensions 

To demonstrate the procedure used for finding the Casimir free energy and 
force we look at the well known case: A massless, scalar field in fiat spacetime 
(Minkowski metric) and no extra spatial dimensions. An additional motiva- 
tion for including this section is that F{Mj~^ = 0) = FMmki as mentioned 
earlier. 

Consider first the general formalism. With DD or NN boundary condi- 

71 

aT 
(2^2^ 



tions we obtain, when making use of the substitution z = a;a-\/M^~+~ef 

t:^DD,NN ^'^ -p/ \ • 

^ = ~ /o \9 r(g) smyrg 



We expand the denominator and use the relation 



_ i).-ig-/.xda; = ^ f-) 'r(z.)K,„.(/i), (49) 



2 
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In the limit s — > — 1 we use the property r(x) sinvrx = 7r/r(l — x) to get the 
free energy for arbitrary T 



:^DD,NN 



^""i:Ei:(^)'-i(-«v^).(so, 

— , Af,, Tl=1 \ / \ / 



(2vr: 



2 



The same expression follows if one makes use of zeta regular izat ion. The 
Abel-Plana formula is powerful in the present context, as it is easily ad- 
justable to different choices for the boundary conditions. 

Consider next fiat space. With Mjv = we obtain from Eq. (ISUj) 

^-•r-if (^)'^i")- 

The first term corresponds to / = 0, and is derivable for instance by taking 
into account the properties of K^{z) for small arguments. 

For high temperatures, aT ^ 1, the expression ( l5Ti) is suitable. The first 
term is the dominant one, as the terms decrease for increasing tempera- 
tures. 

For low temperatures, aT <^ 1, some rewriting is however necessary. We 
go back to the complex function 

rp OO 

F{s) = --ris) Y.^eUkl)-% (52) 

l=—oo kx 

which corresponds to Eq. (139|) when Mjy = 0, s = — 1. Splitting off the / = 
term and using again the Mellin transform fl32l) we can write F{s) as 

rp /•OO 

F(s) = Fi=o-—Y / dtf~^S2Un^TH)e~''''. (53) 

Here 

Fi=o = r(s) $^(fc^)-^ (54) 

and 5*2 (t) is the function 

OO 

S2{t) = Y.e-'\ (55) 

m=l 
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possessing the property [12] 

The first of the three terms coming from the rhs of Eq. (1561) cancels -F)=o, 
leaving 

-i^EE/ *<'-"^'-exp(-«*^-^). (67) 

The first term here is recognized as the zero-temperature energy, F{T = 0) = 
E. With s = — 1 we find 

Furnk = EM^nk ' ^(2^)' Yl ^ {y') ' ^^^^ 



We can now make use of the Abel-Plana formula (H3i) . choosing for the func- 
tion f{z) the form 



fiz) = -(-)'K. (-] 1-y 



(59) 



This leads to, when omitting terms not contributing to the piston model, 

„ OO ^ 

la sr-^ 1 /„„x 

Fuink = EMink - 2^ (4a2^2 ^^2/2^2)2- ^"^^ 

We once more use the Mellin transform, but this time choosing S'2(4a^t) 
together with Eq. (1561) . Some calculation leads to the final expression 

r-,DD,NN t:,DD.NN ^ (^\^ ( 
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where Ej^^^f^ = — 7r^/(1440a^). The corresponding expression for the pres- 
sure is 



pDD,NN _pDD,NN,rj. _ r.^ _ ^ /n \ I / Tlln\ 



n,i=l 

00 



T^\/T/2 v-^ ^2 / T[ln 



a2 



n,i=l ^ ^ 



2 V aT 



(62) 



where Pui^k'^iT = 0) = — 7r^/(480a''). The Casimir energy and force are 
equal to the zero temperature expressions plus correction terms, the latter 
decaying exponentially as T — 0. 

In Eq. ( !62|) we may insert the asymptotic expansion for large arguments, 
K^{z) = (7r/22)i/2g-^[i _^ ^4^2 _ lysz]. It is of interest to extract the 
dominant term in the correction, corresponding to n = / = 1. Approximately 
we then get 

TjDD,NN r>DD,NN/rr n\ , f ^ \ I RO\ 

PKHnk =PMink = ^ ^ """^^ \~ ^) ' ^^^^ 

The physically important point here is that the finite temperature term is 
positive, corresponding to a repulsive force correction (recall that we are con- 
sidering aT ^1). The situation is in some sense similar to that encountered 
in earlier studies when calculating the Casimir force between two parallel 
metallic slabs in physical space, assuming the Drude dispersion relation for 
the material: also in that case the finite temperature effect was found to 
weaken the attractive T = force [H] . 

We now consider the third class of BC's mentioned above: assuming 
Dirichlet boundary conditions on one plate and Neumann on the other we 
find 

The steps are similar to those of the DD and NN calculations, only with 
a factor (—1)" due to the positive sign in the denominator and accordingly 
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^Mink — ~'^/^^Mink ■ "^^^ ^^^^ energy density with DN boundary condi- 
tions becomes 

^ ^ l=-oo Mjv n=l \ / \ / 

(65) 

With Mat = we get 

which is a convenient form for the case of high temperatures. 

For low temperatures we obtain by a similar reasoning as that given 
above, 

F^Zk = E-l, - ^ E (-1)" (y) ^ K. (67) 
(2a)2 ^^^^^ V// ^ \ al J 

with corresponding force 

<i ^p^KT ^ 0) - igj f: (-1)- (^) ' A , (^) 

r'"' (68) 

Again extracting the dominant term by including only n = / = 1 we get 
approximately 

PEZu = PEiUT = 0) - ^ exp (-^) . (69) 

The correction term is the same as in Eq. ( l63l) . but with the opposite sign. 
The thermal correction is attractive. 

5 DD, NN and DN boundary conditions in 
RSI 

Consider first the high-temperature regime. Whereas in flat space this cor- 
responds to aT 3> 1, in RSI the natural choice for high temperatures is 
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T ^ ke'^'^'^". Recall that the lowest values of are ~ ke'^'^'^"; this implies 
aT ^ 1 since ake~^'"^'' ^ 1 for all relevant distances in physical space. In 
this limit Eq. fl50l) is a suitable expression for the free energy and the Casimir 
force is 



^ ' l=-oo Mm n=l ^ 



4 



w E E E '<i { ^"V^'ft + (2-^')- 

^ ^ /=-oo Mjv n=l ^ 



(70) 



After some rewriting we find this in accordance with Eq. (23) in [20]. We 
need only to include the E{M]^ = 0) term to get the Casimir force for a 
massless scalar instead of a massive. 

We find that the high temperature limit is valid for T ^ lO^^K. Only 
temperatures much less than these are expected to be of physical importance. 
It is most natural therefore to find the Casimir energy and force for T <C 
^g-fc7rrc_ Note that the brane low-temperature condition does not fix the 
magnitude of the product aT relative to unity. With k ~ 10^^ GeV, e'^"^"^" ~ 
10~^^ we only get the weak condition T <^ 10^ GeV. As an example, choose 
T = 300 K (2.6 X 10"^^ GeV), a = l/im, from which it follows that aT = 0.15. 
In most cases of practical interest we will have aT -C 1, although one can 
easily consider cases where aT ^ 1, still compatible with the condition 
T < ke-''^^^. 

Using the same procedure as in flat space we find the RSI equivalent to 
Eq. dSH]), 

F... = E,s, - j:j:± (s±i«) ' (iy^) . 

(71) 

We can differentiate this expression to find the Casimir force. By assuming 
dkx/da = —kx/a we get Eq. (17) in [21], only missing the first term. The 
assumption holds for all k^ proportional to 1/a, which is the case for DD, 
NN and DN BC. However, we are using the piston model, and a term like 
the first one of Eq. (17) in [21j should not occur in the Casimir force. The 
force term in question is independent of a, and thus corresponds to a free 
energy term that is linear in a. As a consequence of using the piston model. 
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all terms independent or linear in a must be removed from the free energy. 
Although the metric in [2T] does not include the warp factor e~'^^^'''^ the 
expressions are the same, since the warp factor only affects the values of the 
Mivs. In Eq. flTTl) Ersi is the zero temperature energy in RSI and is found 
from Eq. (HOj) using the Abel-Plana formula (H3|l with 



After some variable changes the energy reads 

and for DD and NN boundary conditions it simplifies to 

= -si; £ E (74) 

ra=l Mjv 

The Casimir force at zero temperature is 



poo,™,^ = 0) = - ^ 5: 5: (2aM„„) 

n=l Afjv 

n=l Afjv 



(75) 



This is in accordance with [2B] through that paper does not consider the 
Casimir effect rising from a bulk scalar in the RS model in particular. Insert- 
ing the approximation Eq. (128|1 for M^v we see that the energy is essentially 
the same as in [10] . There are three minor differences. First of all the energy 
in [in] has some extra terms linear and independent of a since the piston 
model is not used. Secondly, factors p are included to make the expression 
hold for electromagnetic fields, where p is the polarizations of the photon. 
The last difference is a factor of 2 included to account for 'the volume of 
the orbifold'. Since we can not see how this factor occurs, it is not included. 
This is also equal to Eq. (26) in [20]; only it contains the E{M]y = 0) term 
since a massless field is considered. 
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The summation over in Eq. fl7ip is still left and can be done using the 
Abel-Plana formula with 



f{z) = - (M^ + z'/a') if 3 i-^zya^ + M\ 



TT 



After inserting Kz^z) = (7r/22;)^/^e ^(1 + 1/z) the free energy reads 



(76) 



""rsi — Ersi 



EE'-' 

Mjv 1=1 
l_ 

T 



dz- 



^j=0,a l-{l3jz/a)2 



(ft)/a-l)(/3./a-l) 2z 
(A,/a+l){/3„/a+l)^ 



1 



sm 



T 



(^/a)2 - 



(77) 



We continue by inserting the /3's for DD and NN boundary conditions, ex- 
pansion of the denominator and the variable exchange x = z/a to get 



Th 



^RSI — ^RSI ^2 



EE'-'E 

Mn 1=1 n=l 



dxe 



-2nax 



T 



— cos 



T 



I 



x2 - Ml 



(78) 



Integrals of this form is solved in Appendix A with the result 

I dx (^sin (^AVa;2 - C^^ - A^x^ - cos [a^Jx^-C^^) ^'^"^ 

giving the free energy 



(79) 



pDD,NN _ pDD,NN ^ \ ^ \ ^ \ ^ 

-^flS/ -^RSI ^2 

Mn 1=1 n=l 



Ml 



Z-.^^(2nay + (l/Ty 



(80) 



X K2 (MjvV(2na)2 + (//r)2 
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and force 
pS^"" = PSr^'iT = 0) 



-M^ttti ii2nay + il/m' ^ (81) 

-AEEE ""^^'^^^^ 3 i^.(MW (2na)2 + (VT)2 ). 
^ t^tr^((2na)2 + (//T)2)5 V y 

This expression, belonging to the low-temperature regime T -C ke~^'^^'', can 
be used both for aT ^ 1 and for aT ^ 1. The argument of the Bessel 
functions will always be large since ake~^'^^'' ^ 1 for all relevant distances. 
The correction terms to the zero temperature energy and force expressions 
are small. The expression has to our knowledge not been given before. The 
leading term for the force in terms of T/M^v is 

pDD,NN pDD,NN/ 




-./ (2anT)2+P / T ^ (2naT)2 (82) 



T 



In contrast to fiat space, both zero temperature Casimir force and the thermal 
correction is negative. Hence the Casimir effect in RSI is stronger in the low 
temperature limit (T <^ ke''^'^^") both for aT -C 1 and aT ^ 1. 

The DN expressions deviate from the DD and NN expressions in RSI 
in the same way as in flat space. The factor of (—1)"" must be included 
in sum over n, where the sum over n originates from the expansion of the 
denominator in Eq. (173!) and (!77|) . 



5.1 Comparison to flat space 

The reason for calculation the Casimir force in RSI is to find out where there 
are deviations from the Casimir force in fiat spacetime without extra spatial 
dimensions. For an easier comparison we give the full expression for the 
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Casimir force of a massless hulk scalar m RSI with DD/NN BCs. 



P. 



DD,NN 
RSI 



TT 



480a4 



^ oo 



3 (aM. 



N 



-K2 {2aMNn) + 



1 (aM, 



N 



471^ n 



ra,i=l 



3(aT)2 / vr/n \ /— — —74^ 



^ 00 



M]v ra,;=l 



v/2 v/y "3 

(aMjv)=^(2n)^ 
((2n)2 + (//aT)2)i 



n2 / Tr/n 

5 



7/ 

i^i faMjvV(2n)2 + (//Ta)2 



(83) 



7^2 faM,v\/(2n)2 + (//aT)2 



(aM^)2 (3(2n)2 - (//aT)- 
((2n)2 + (//aT)2)2 



This expression is good at low temperatures because the Bessel function 
decreases exponentially at high arguments. Hence, we only need to sum over 
the first couple of values from M^v, n and / if the other factors (M^r, a and 
T) ensure that the argument is much greater than one. On the other hand, 
if the argument of the Bessel function is not large we need to be careful that 
the sum has converged. Now the essential question is: For what values (of 
Mtv, a and T) are the sums of Bessel functions are of the same magnitude 
as the flat spacetime at zero temperature (i.e. the flrst term)? Or simply: 
When can we see a deviation from the ordinary Casimir force? 

Looking at Eq. (15^ we see that at zero temperature we need aMf^ ~ 1 
for noticeable difference. We know that Mm ~ ke~^'^^'' for low N, and k is 
usually set to ~ Mpi ^ lO^^GeV in RSI. In the original paper of Randall and 
Sundrum they propose to choose krc ~ 10 in order to solve the hierarchy 
problem. With these values we flnd that a is ~ lO^^^m. There is no point 
in looking at distances smaller than the size of an atom. Only distances of 
physical relevance (> Inm) are of interest. In Figure[2]we keep k = lO^^GeV, 
but choose e~^'^^'= = lO"^^. xhe difference from RSI to ordinary Casimir force 
FMink at the zero temperature is given in Fig. [2j By choise of parameters 
the magnitude of the correction in RSI is of the same order of magnitude 
as FMink, given by the red line in Fig. [2] (a). With smaller value of krc we 
will not se any difference at separations larger than Inm. The corresponding 
size of the extra dimension is rc ~ 10~^^m. In Fig. [2] (b) we se the ratio of 
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i-igure A. i^Rgj 



-P^Hnk^ for a massless scalar bulk field with k 



lO^^GeV 

and d = e"'^'^^'' = 10^^^. Left panel: The difference between the Casimir 
force for RSI and fiat 3+1-spacetime at T = and T = 3 x 10^ K. The 
Casimir force for a massless scalar field in flat spacetime at zero temperature 
is included for comparison. Right panel: The ratio of the difference between 
the Casimir force in RSI and fiat spacetime at T = and T = 3 x 10^ K to 
the force at zero temperature in fiat spacetime. 



this difference to the Casimir force at zero temperature. We observe that the 
extra term we get in RSI, goes faster to zero than Fuink- Now we turn to 
relevant values of the temperature. The choices of rc, k and a are still the 
same, and after some testing it turns out that a^axT ~ 3 is suitable. To be 
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sure that the sums have converged we let both n and / run to 30. The result 
is presented in Fig. [21 The green line is the Prsi — Puink for T = 3 x lO^K 
and we see that this gives a stronger Casimir force than at zero temperature. 



5.2 Comparison to flat space with one extra dimension 

In the previous section we compared the Casimir force in RSI with the 
Casimir force for a massless scalar i ordinary fiat 3+1-spacetime. In this 
section we will look at a higher dimensional spacetime which is fiat, i.e. no 
warp factor in the metric. As mentioned in the introduction, this topic has 
gained a lot of interest lately. The extra dimension is a torus with circum- 
fence 2tiL. In this case = N/L, with = 0, ±1, ±2, . . .. With these 
values for the Mats instead of the ones we have in RSI we see that Eq. fITOl) is 
equal to the high temperature expression in |19] . However we have not found 
an expression corresponding to Casimir force in Eq. fl83|) . Hence, we will 
derive such an expression using the Chowla-Selberg formula in [45] . With 
the new values of the M/vs, the function F{s) (for DD BCs) reads 

~ oo oo oo 

F{s) = --V{s)Y, E E((2vrTr+(W«)' + (W)-^ (84) 

l=—oo N=—oo n=l 

Rewriting this to homogenous Epstein zeta functions 

oo ' 

ZE,p{s;ai, . . . ,ap) = ^ {{aikif + {a2k2f + . . . + {apkpf)~\ (85) 

fci,...,fcp=— oo 

we find 

F{s) = -I^V{s)Ze,2{s; 27rT, l/L) - I-V{s)Ze,^{s; 27rT, vr/a, 1/L). (86) 

The notation A:p=-oo Cleans that for ki to kp we sum from — oo to oo. 

The prime / behind the sum means that the term k\ = k2 = . . . = kp = Q is 
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omitted. The Chowla-Selberg formula is 



TT^r (s - f ) 



+ 



+ 



ZE,m{s; ai, . 

m 



(n: 



ai 



r s) 



■' E,p—m I S 



2 ) Q'm+li 



r(3) (n: 



^1 - -- -tAjT^ 



(87) 



-m/2 



\ i=i \ 



p 

j=r?i+l 



After use of this with m = 2, we put s = — 1 and remove all terms linear and 
independent of a. Then we see that the free energy is 



T 



T{-1)Ze,2{- 

lOTT 

oo ' oo 
d ^ — V — > 



N= 



-oo n,Z= 



1; 2ttT, it /a) 
{2naf + {l/Ty 



K2 [N/L^{2naY + {l/TY 



(88) 



The first term can be identified as F{Mj\f = 0) and the second term is equal 
to Eq. (IHOl) (when Mjy = N/L). From Yl^i=-oo factor 4 when 

we let the sums over n and I run from 1 to 00. The term / = 0, but n 7^, 
gives the zero temperature expression and n = with Z 7^ is independent 
of a and should be removed in the piston model. 

Since Eq. f lHOl) leads to Eq. fl83|) we can conclude that we get the same 
answer with the Chowla-Selberg formula as with the Abel-Plana formula in 
flat spacetime with one extra spatial dimension. However, the Abel-Plana 
formula can be used regardless of the values of M^v, while the Chowla-Selberg 
formula is only useful when we can rewrite our expressions to homogenous 
Epstein zeta functions. The second advantage of the Abel-Plana formula is 
that different boundary condititions can easily be obtained. 



25 



6 DD, NN and DN boundary conditions in 
RSII 

In RSII the Kaluza-Klein modes are continuous and we must replace the sum 
over Mtv with an integral, 

In Eq. (ISOll we get an integral on the from 



dz {z^ + A^) ^ Ks [b^z^ + = ^J^A^K2{AB). (90) 

The derivation of this formula is given in Appendix A. We find that the free 
energy and force in RSII are 

rp 3 rT-i3 OO OO ,2 

i.DD.NN ^ ^ _^yy L^K, (AnTaln) (91) 

^^^^ lUOka^ ak ^ ^ ^ ' ^ ' 

1=1 n=l 



and 



^^^^ lUOka^ a?k ^ ' 

[=\ n=l .g2) 

l=\ n=l 

These expressions are convenient for the high temperature limit aT ^ 1. In 
RSII there are only two temperature regimes, aT ^ 1 and aT -C 1, since 
Mtv is continuous. To find the low temperature limit we insert Eq. ( I89ll into 
Eq. flHOj) and use the integral 

d^ z'KM^) = ^. (93) 

The free energy reads 

FnSir = Ei^iT -^yy r, (94) 

2vrfc^^((2an)2 + (//T)^)t' 
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where E^^j^^ = — j^^^ is the zero temperature energy in RSII. We can 
find this energy from e.g. Eq. (|7il) by making use of Eq. (l93ll . We use the 
Melhn transform as in low temperature, fiat spacetime (with 6*2(40^^)) and 
find ^ 

Fnsn = Ensn ' ^2 1^ [j) ^2 ■ (95) 

n,l=l ^ ' 

The Casimir force is 

Pnsn =Pnsn (T = 0) + -K, ^— , (96) 

n,l=l ^ ^ 

with pDD.NN.^ - - 3Cfi(5) 

In the low temperature limit aT <^ 1 we get from Eq. fl96l) the dominant 
term corresponding to n = / = 1, 

p-,r^«r(r^O).^;fexp(-^). (97) 

The temperature correction term is repulsive. 

As in RSI the only difference between DD/NN BC and DN is a factor 
(— 1)" in the sum over n. 



7 Concluding remarks 

Our main objective has been to calculate the finite temperature Casimir effect 
for a scalar field residing in the bulk in the two Randall- Sundrum models, 
RSI and RSII. Two parallel plates are envisaged, with gap a, located on one 
of the RS branes. We have given most attention to the RSI model. Robin 
boundary conditions, cf. Eqs. (fT3|l and ([HD, are assumed on the two plates. 
The geometrical picture is the piston model, as illustrated in Fig. 1. We have 
made use of the Abel-plana summation formula throughout, as this turns out 
to be the most convenient choice in the present context. 

In the case of fiat space the basic expressions for Casimir free energy and 
force (per unit surface area) are worked out in the form of series in Sect. 4, 
both for high and for low temperatures. A characteristic feature for DD and 
NN boundary conditions on the two plates is that the dominant part of the fi- 
nite temperature correction term for low temperatures {aT ^ 1) is repulsive. 
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That is, the force becomes decreased shghtly when the temperature increases 
from zero. In this sense the behavior is analogous to that encountered in the 
case of conventional Casimir theory for metallic slabs in physical space when 
the dispersive relation for the material is taken to have the Drude form |31] . 

The RSI model is covered in Section 5 in an analogous way. The dominant 
term in the Casimir force shows the characteristic property to strengthen 
the zero temperature effect instead of weaken it as in RSII. From Eg. (1811) 
we can evaluate the Casimir force both for aT <^ 1 and aT ^ 1, provided 
T <^ ke~^'^^''. This section also covers a comparison to flat space with and 
without a compactified extra dimension. 

In Section 6 the RSII model is considered. We have dM/ k = Tr/k dM/ {2tt) 
and thus the Casimir force has a characteristic n/k times the Casimir force 
of a 4+1-dimensional Minkowski spacetime. This is pointed out by Morales- 
Tecotl et al. |T2], but as an argument against using zeta functions in the 
regularization and rather use Green's functions. It is not only the regular- 
ization method that is different, the physical picture is also differs from this 
article and the work by Frank et al. [TOl [H] . While Frank et al. calculate 
the free energy of a slice of the bulk, Morales-Tecotl et al. try to restrict the 
system to the brane by evaluating the Green's function at the visible brane 
(?/ = for RSII). In this way they claim to incorporate the localization prop- 
erties of the modes of the scalar bulk field (ipNiu))- Note that the Green's 
function method includes a integral over y and by setting y = 0, Morales- 
Tecotl et al. thus remove the y-dependence of a part of the integrand before 
integrating. However, keeping the y-dependence before the integration makes 
one obtain results different from those of Frank et al. The delicate point is 
how to include the localization properties of the modes. In Ref. |36] it is 
proposed to resolve the issue by changing the boundary conditions. The 
problem with the localization properties of the modes needs to be resolved 
before the Casimir effect form an electromagnetic field can be considered. 
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A Integrals and Bessel functions 

We need to calculate the nontrivial integral 

f dx (^AVx^ - C2 cos (^AVx2 - - sin (^AVx^ - ^2)) e'^^ (A.l) 



First use the substitution u = \/ x'^ — C"^ to find 

" -^^={Au^ cos(Am) - u sin(An))e-^^"'+^' 

=A^\ du ^ ^ cos(A^)e-^^^ 

/•oo 1 

-A / dw , cos(AM)e-^^^^^ 

io Vu2 + C2 



(A.2) 



du 



+ C2 

Using Eq. 3.961 in 



u sin(y4u))e" 



/•°° ^^g-/3v/^sin(ax) = ^ i^i (A.3) 

as well as 



'0 

we get 



/" da; (^AVx2 - C2 cos (^AVa;2 - C2 j - sin (^AVa;2 - ^2) j e~^^ 

- AM^Kq {C\/A^ + 52) . 

After differentiating and using the relationships between K^^ K\ and K2 we 
find the solution 



I dx (^Av/;^?2r^cos (^AVx2 - C2 j - sin (^AVx2 -C2jj 
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Next consider the integral 

' (A.7) 



With the substitution m = a/ (-zM)^ + 1 we get 



— (A.8) 



25 y^I^ V 

again using Kziyz) = (tt/2zY^'^ e^^{l + l/z). As the integral representation 
ofK2{BA) is I?]) 

i / 1 \ 2 

^2(ax) = ^^^^!^ du e— - (A.9) 
we find after several partial integrations 

Thus the integral becomes 

dz{z^ + A^f^Kz [b^Jz^ + = ^J^A^K2{BA). (A.ll) 
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